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Specific lieat and ac magnetic susceptibility measurements, spanning low temperatures (T > 
40 mK) and high magnetic fields {B < 14 T), have been performed on a two-dimensional (2D) anti- 
ferromagnet Cu(tn)Cl2 (tn = C3H10N2). The compound represents an S = 1/2 spatially anisotropic 
triangular magnet realized by a square lattice with nearest-neighbor {J/ks = 3 K), frustrating next- 
nearest-neighbor (0 < J'/ J < 0.6), and interlayer {\J"/J\ ~ 10"'') interactions. The absence of 
long-range magnetic order down to T = 60 mK in B = and the behavior of the specific heat 
for T < 0.4 K and B > are considered evidence of high degree of 2D magnetic order. In fields 
lower than the saturation field, Bsa.t = 6.6 T, a specific heat anomaly, appearing near 0.8 K, is 
ascribed to bound vortex-antivortex pairs stabilized by the applied magnetic field. The resulting 
magnetic phase diagram is remarkably consistent with the one predicted for the ideal square lattice, 
except that Bsa,t is shifted to values lower than expected. Potential explanations for this observation, 
as well as the possibility of a Berezinski-Kosterlitz-Thouless (BKT) phase transition in a spatially 
anisotropic triangular magnet with the Neel ground state, are discussed. 

PACS numbers: 75.40.-s, 75.10.Jm 



I. INTRODUCTION 

Two-dimensional quantum antiferromagnets have at- 
tracted a significant amount of theoretical and exper- 
imental attention due to the unconventional magnetic 
properties resulting from the interplay between quantum 
fluctuations and geometrical frustration.^"'^ One example 
is the S = 1/2 spatially anisotropic triangular antiferro- 
magnet, which can be treated as a square lattice with the 
nearest-neighbor (nn) interaction J and frustrating next- 
nearest-neighbor (nnn) interaction J', Fig. 1. Between 
the limiting values of J' = (ideal square lattice) and 
J' / J 3> 1 (spin chain), several phases exist for varying 
ratios of J'/J,"*"^ while the presence of an applied mag- 
netic field provides an additional constraint. Consider- 
able theoretical interest in the spin liquid phase [J' / J > 
1) in a magnetic field^"^^ was triggered by experimental 
studies of CS2CUCI4, a spatially anisotropic triangular 
magnet with J' /J w S.^^-w 

Recently, Cu(tn)Cl2 has been identified as a poten- 
tial model system for the realization of the spatially 
anisotropic triangular lattice from the coUinear Neel 
phase (J'/J < 0.6)}^ For Cu(tn)Cl2 studied in B = 0, 
no evidence for long-range magnetic order was observed 
down to 60 mK, and the data suggested intralayer inter- 



action strengths of J/ks = 3 K and < J'/J < 0.6, 
while the interlayer coupling is \ J" / J\ « 10"'^. These 
interactions are described by the Hamiltonian 

■H = J^S,;-S, + J'^S,-S,+ J"^S,-Sfc, (1) 

nn nnn i,k 

where i,j label intralayer spins and k labels the interlayer 
ones. 
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FIG. 1: Realization of Heisenberg model of a spatially 
anisotropic triangular lattice within a single be layer in 
Cu(tn)Cl2. The layers are stacked along the a direction. The 
full circles denote Cu^"*" ions. 
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The motivation of the present work was to explore the 
response of Cu(tn)Cl2 in _B ^ 0, especially at low tem- 
peratures, T <C J/Ub- For this purpose, experimental 
specific heat and ac susceptibility studies were performed 
over a wide range of temperatures (40 mK < T < 10 K) 
and magnetic fields (0 < S < 14 T). The low value 
of the intralaycr exchange coupling affords easy access 
to the magnetic phases below and above the saturation 
field. On the basis of the magnetic field induced fea- 
tures observed in the specific heat and ac susceptibility, 
a magnetic phase diagram is constructed and analyzed 
within a model of a field induced Berezinski-Kosterlitz- 
Thouless (BKT) phase transition theoretically predicted 
for the pure 2D square lattice with zero nnn coupling. 
It is noteworthy that the long-range orderings in some 
quasi-2D, square lattice systems, namely Sr2Cu02Cl2 
and several Cu(pyz)2-based materials (pyz = pyrazine), 
have been placed close to the BKT transition. ^^"^^ 

Our presentation begins with a discussion of the sam- 
ple and experimental details, which are followed by a 
description of our experimental results. Next, an anal- 
ysis and a discussion of the results are given before the 
salient points are assembled into the magnetic phase di- 
agram. The paper concludes with a summary and some 
comments about possible future directions. 



II. SAMPLE AND EXPERIMENTAL DETAILS 

The crystal structure of Cu(tn)Cl2 (tn = C3H10N2), es- 
tablished at 150 K, is orthorhombic (space group Pna2i) 
with the lattice parameters a = 17.956 A, 5 = 6.859 A, 
and c = 5.710 A.^^ The structure consists of covalently 
bonded ladders running along the c-axis, while the ad- 
jacent ladders in the 6c-planc are linked through inter- 
molecular N— H- • • CI hydrogen bonds formed by all four 
H atoms of the amino groups. In the a direction, the 
layers are connected by weak C— H- • • CI type interac- 
tions. The strongly elongated octahedra coordinating 
the Cu(II) ions stabilize the d^i electronic ground state. 
Consequently, the propagation of exchange pathways be- 
tween Cu(II) ions leads to the formation of a spatially 
anisotropic triangular lattice in the &c-plane (Fig. 1). 

The synthesis of all Cu(tn)Cl2 samples followed the es- 
tablished procedure, which produces polycrystals that 
are powdered and pressed into pellets (nominally 3 mg to 
10 mg) for the specific heat studies or placed into appro- 
priate specimen holders for the magnetic susceptibility 
investigations. A sample with a mass of 26 mg was used 
in the ac magnetic measurements. 

Using several experimental probes and instruments, 
the specific heat measurements were performed over the 
temperature range from 100 mK to 10 K and in magnetic 
fields up to 14 T. More specifically, the studies in the 
millikelvin temperature region and in magnetic fields up 
to 2.5 T were performed using a relaxation calorimeter^" 
mounted on a dilution refrigerator. A commercial (Quan- 
tum Design PPMS) device was used for the specific heat 



studies over the temperature range from 1.8 K to 10 K 
and in fields up to 9 T, while the low temperature stud- 
ies down to 0.35 K and in fields up to 14 T utilized an- 
other commercial instrument equipped with a ^He insert. 
In each instance, the contribution of the background ad- 
denda was determined in separate runs. Finally, the sepa- 
rate study of the specific heat of a diamagnetic isomorph, 
Zn(tn)Cl2, allowed the phonon contribution to be deter- 
mined. 

The magnetic susceptibility studies were performed 
with ac (232 Hz) mutual inductance coils mounted on 
a dilution refrigerator equipped with a 10 T magnet. 
With the sample immersed in pure ^He that provided in- 
timate thermal contact with the mixing chamber, the in- 
phase and out-of-phase signals of the susceptibility were 
recorded by a two channel lock-in amplifier. Typically 
the data were obtained by isothermal field sweeps at a 
rate of 50 mT/min, and the data were independent of the 
direction of the field sweep. 



III. EXPERIMENTAL RESULTS 

A. Specific heat in S ^ 

Evidence of the presence of 2D magnetic correlations 
was initially observed as a round maximum in the tem- 
perature dependence of the total specific heat, Ctot(7'), 
near 2 K in i3 = O,'^^ and our present work found this fea- 
ture to evolve with magnetic fields up to 9 T, Fig. 2. In 
fact, the influence of the magnetic field can be separated 
into two regimes. For _B < 4 T, the height of the spe- 
cific heat maximum, Cmaxj decreases with increasing field 




T(K) 



FIG. 2: Temperature dependence of the total specific heat, 
Ctot, of Cu(tri)Cl2 in various magnetic fields. Solid lines 
represent data for B = G, 0.5, 0.75, 1, 1.5, 2, 2.5, 3, 3.5 and 
4 T, while the dashed lines correspond to -B = 5, 6, 7, 8 and 
9 T. The dot-dashed line represents the specific heat of the 
diamagnetic isomorph Zn(tn)Cl2. 
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while its temperature, Tmax, remains nearly unchanged, 
whereas for B > 4 T, Cmax increases with increasing 
field and Tmax shifts towards higher temperatures. The 
phonon contribution to the specific heat can be approx- 
imated by the specific heat of the diamagnetic isomorph 
Zn(tn)Cl2, Fig. 2. These high temperature measure- 
ments were extended to lower temperatures in i? < 2.5 T, 
and the results of the total specific heat are shown in 
Fig. 3. These studies revealed the presence of an anomaly 
appearing at about 0.8 K when B ^ 0, and this feature 
develops with increasing magnetic field (inset of Fig. 3). 
The magnetic field and temperature dependences of this 
anomaly were measured in i? < 14 T, Fig. 4, where this 
feature reaches a maximum in 4 T at 0.8 K and decreases 
in magnitude, with a simultaneous shift to lower temper- 
atures, until its presence is no longer resolved in fields 
greater than 7 T. 




T(K) 



FIG. 4: (Color online) Temperature dependence of the total 
specific heat of Cu(tn)Cl2 in various magnetic fields. 



B. Magnetic susceptibility in i? < 10 T 

The results of isothermal ac susceptibility studies in 
i3 < 10 T are shown in Fig. 5, and no significant hystere- 
sis was observed between the up and down sweeps. The 
data do not possess any sharp anomalies corresponding to 
a phase transition, but the appearance of a shoulder near 
6 T corresponds to the field where the low temperature 
specific heat anomaly vanishes. Unlike the low temper- 
ature specific heat anomaly, the shoulder survives up to 
1 K, suggesting that it is associated with the saturation 
magnetic field i?sat- 
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FIG. 3: Temperature dependence of the total specific heat 
of Cu(tn)Cl2 in -B = 0, 0.5, 0.75, 1, 1.5, 2 and 2.5 T. Inset: 
Temperature dependence of the difi'erence AC = [C{T,B) — 
C{T, 0)] between the total specific heat in finite and zero mag- 
netic fields. 



FIG. 5: (Color online) Magnetic field dependence of isother- 
mal ac susceptibility of Cu(tn)Cl2 at various temperatures. 
Inset: Magnetic field dependence of the normalized isother- 
mal magnetization at 40 mK. 



IV. ANALYSIS AND DISCUSSION 

A. Magnetic correlations in i? = 

As stated earlier, previous specific heat studies in 
B ~ did not indicate any magnetic phase transi- 
tion down to 60 niK, and a clear quadratic dependence 
was observed at low temperatures. The latter coincides 
with the expected 2D character of short-range magnetic 
correlations,^^ since the spin wave analysis of low dimen- 
sional models with nn interactions predicts a behav- 
ior for the low temperature specific heat of square and 
triangular lattices and a T dependence for a linear chain 
model. ^"'^ The absence of a A-like anomaly associated with 
magnetic ordering is another intriguing feature of the spe- 
cific heat data. 
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In the absence of the frustration (J' = 0), the order- 
ing temperature Tn of the isotropic square lattice with 
interlayer coupUng J" can be expressed as^^'^^ 



(2) 



where M/Mq is a staggered magnetization, a represents a 
lattic:c constant, and ^ is an intralayer correlation length 
given by 



a 



= 0.5 



exp 



{ksT ) 



1 



27rp,, 



(3) 



Here, ps is a spin stiffness, and for the isotropic square 
lattice with intralayer nn exchange coupling J, it can 
be expressed as ps = 0.18J. Using the parameters 
\J" /ks] ~ 3 mK, iJ/fcsl 3 K, determined in Ref. 15, 
and (M/Mq)^ « 0.3, a phase transition in Cu(tn)Cl2 
might be expected at Tn » 0.8 K. The absence of the 
phase transition down to 60 mK may suggest a signif- 
icant reduction of both the spin stiffness and the stag- 
gered magnetization as a consequence of frustrating nnn 
J' coupling.^"^'^^'^^ Apart from the aforementioned re- 
duction of staggered magnetization, the weakness of the 
interlayer coupling can also lead to the absence of a de- 
tectable phase transition at finite temperatures in spe- 
cific heat measurements. Recent Monte Carlo studies of 
the specific heat of an isotropic square lattice with vari- 
ous interlayer couplings revealed that, for J" / J < 0.015, 
the peak associated with the phase transition vanishes.^'' 
Since \J"/J\ « 10~^ for Cu(tn)Cl2, the ordering effects 
might not be observed in the specific heat even in the 
absence of frustration. While experimental studies of 
Cu(pyz)2 (0104)2 are consistent with this scenario, the 
observation of a phase transition in the specific heat of 
Ou(H20)2(02H8N2)S04 with a comparable J"/ J ratio 
shows the need to consider additional effects. 




FIG. 6: Magnetic field dependence of the value for the maxi- 
mum of the Cu(tn)Cl2 magnetic specific heat, Cmax, divided 
by the gas constant R (open squares). The solid line is a guide 
for eyes. Magnetic field dependence of the temperature, Tmax, 
denoting the position of Cmax (full squares). The broken line 
is a linear fit (see text). 



Using g = 2.12 and J/ks = 4 K, which correspond to the 

best estimates of the parameters resulting from the lin- 
ear chain model as applied to 0u(tn)0l2 in Ref. 15, Eq. 4 
predicts Bgat = 5.6 T, a value that coincides rather well 
with the crossover region appearing between 4 T and 6 T. 
However, in comparison with the linear chain, which is 
characterized by a decrease of Tmax with increasing field 
for fields lower than Sg^t and an increase for higher fields, 
the observed shift of Tmax with respect magnetic field be- 
haves differently. Alternatively, the effect of an external 
magnetic field on the short-range correlations and the 
thermodynamic properties of a square lattice have been 
theoretically investigated,^^ and the saturation field is 



9MB.Bsat = 4 J 



(5) 



B. Magnetic correlations in B 7^ 

1. Specific heat for T > J/ks 

After the subtraction of the phonon contribution, ap- 
proximated by the specific heat of the diamagnetic iso- 
morph Zn(tn)Cl2 (Fig. 2), two different regimes can be 
distinguished in the behavior of Cmax for T > 1.8 K. 
For _B < 4 T, Cmax decreases with increasing field while 
Tmax remains nearly unchanged, whereas for i? > 5 T, 
Cmax increases with increasing field and Tmax 

shifts to- 
wards higher temperatures. Fig. 6. Such a qualitative 
behavior has been predicted for a linear (ID) chain with 
nn coupling J",^"^'^^ where the crossover between the two 
regimes occurs at a saturation field 

5MB-Bsat = 2J . (4) 



For 0u(tn)0l2, this model predicts i?sat = 8.4 T, which is 

too high to correspond to the observed crossover region. 

This qualitative comparison of the limiting ID and 
2D theoretical models with the experimental data sug- 
gests the existence of two phases in 0u(tn)0l2, namely 
a low field phase stable in the magnetic fields below 
the crossover region where antiferromagnetic correlations 
prevail and a high field paramagnetic phase stabilized in 
the fields above the crossover region. However, the atyp- 
ical magnetic field dependence of Tmax suggests a more 
complex low field phase. Indeed, a closer look at the 
behavior shown in the Fig. 6 reveals that, within experi- 
mental uncertainty, Tmax remains constant for _B < 2 T, 
before increasing slightly in the fields from 2 T to 4 T, 
and then experiencing a quasi- linear increase for the fields 
above 4 T, as demonstrated by a linear fit performed in 
the range from 4 T to 9 T. 
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FIG. 7: CmT vs. r dependence for Cu(tn)Cl2 in B = 0, 0.5, 
0.75, 1, 1.5, 2 and 2.5 T. Inset: Magnetic field dependence of 
the coefficient 6, Eq. 6. The solid line represents a linear fit. 



2. Specific heat for T < J/ks 

The quantitative analysis of magnetic specific heat, 
Cm, at temperatures T < 0.4 K and B < 2.5 T re- 
veals the dependence, as can be seen in Fig. 7, where 
the data are plotted as CmT^ vs. T^. This approach as- 
sumes the magnetic specific heat can be expressed as a 
sum of a/T^ and bT'^ terms, where the first term corre- 
sponds to the contribution of the nuclear spins and/or 
long-range correlations between electronic spins, and the 
second term is associated with 2D short-range correla- 
tions. A set of linear fits of the individual CmT"^ vs. 
dependences for fixed magnetic field was performed in the 
temperature range from nominally 150 mK to 400 mK. 
The fitting procedure revealed a monotonic rise of the co- 
efficient h (in imits of J/(K'^ mol)) with increasing mag- 
netic field and a « 0. The monotonic increase can be 
approximated by a linear dependence 



b{B) = 1.36 -I- 0.30 B 



(6) 



as shown in the inset of Fig. 7. The observed low tem- 
perature behavior of Cm (T) siiggcsts the preservation 
of the mainly 2D character of the magnetic correlations 
in B ^ 0. Furthermore, the field does not introduce an 
energy gap to the excitation spectrum, which remains 
gapless as expected for the square lattice in magnetic 
fields below the saturation value. 

To more precisely trace the position of the low tem- 
perature peak, Tp, appearing at about 0.8 K, as a func- 
tion of magnetic field, the specific heat in zero magnetic 
field was chosen as a reference background that was sub- 
tracted from the specific heat in nonzero magnetic field, 
inset of Fig. 8. Using the low- field data obtained in the 
same way (inset of Fig. 3), the magnetic field dependence 
of Tp can be extracted (Fig. 8). It can be seen that a 
naive extrapolation to zero field provides Tp « 0.7 K, 
which is close to the value of the transition temperature 



FIG. 8: Magnetic field dependence of the peak position, Tp, 
of the low temperature specific heat anomalies (full squares) 
in Cu(tn)Cl2. The theoretical prediction for the field induced 
BKT phase transition for the isotropic square lattice is de- 
noted by the open squares and the broken line. Inset: Tem- 
perature dependence of the difference between the specific 
heat in finite and zero magnetic field. The solid lines repre- 
sent _B = 1, 2, 3, 3.5 and 4 T, and the dashed lines correspond 
to B = 4.5, 5, 5.5, 6 and 7 T. 



Tn ~ 0.8 K estimated for the isotropic square lattice. 
This coincidence might support the suggestion that the 
weakness of the interlayer coupling prevents the obser- 
vation of the phase transition in the specific heat. An 
extrapolation of the high-field experimental data yields 
S(Tp -> 0) « 6 - 7 T, which is lower than B^^t = 8.4 T 
calculated for the isotropic square lattice. 

Finally, to determine i?sai from the specific heat data, 
we analyzed the CmI?") data in 9 T and 14 T. While 
the peak height of Cm{T) in 9 T is still much lower than 
Cmax/R = 0.438, the theoretical prediction for the S = 
1/2 ideal paramagnet (Fig. 6), the corresponding B = 
14 T value of Caia^/ R = 0.42 is close to the expected 
result. For B > Bgat, a gap A opens in the spin excitation 
spectrum, developing linearly with magnetic field as^^ 



A = 5Mb(S-Ss, 



(7) 



In the temperature region where thermal fluctuations 
overcome the interlayer coupling and T < A/fcs, a 2D 
character of magnon spectra can be expected, resulting 

13 



m 



Cm{T) « - exp(-A/fcBT) 



(8) 



Fitting the Cm(T < 1 K, B = 9 T) and Cm(T < 3 K, 
B = 14 T) data yields A/fcs = 3.4 K and 10.5 K, respec- 
tively, and it follows from Eq. 7 that Sgat = 6.6 ±0.1 T, 
when assuming g = 2.12. It is important to note that the 
significant differences between the observed Sgat = 6.6 T 
and the 8.4 T value expected for the isotropic square 
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lattice indicate that Cu(tn)Cl2 is indeed affected by the 
presence of frustrating nnn interactions, since the inter- 
layer coupling is expected to increcise -Bg^t cis 



5MBBsat = 4J + 2J" 



3. AC susceptibility for T < J/ks 



(9) 



The ac susceptibility data, x{T,B) = dM/dB, can 
be integrated to obtain the magnetization M{T,B), as 
shown in the inset of Fig. 5, and M{T 0, B) data 
is commonly used to establish Bgat, if the saturation 
plateau is accessible. However, in many instances, the 
transition to the fully-polarized state is broadened by 
several effects, including issues related to finite temper- 
ature, orientation of the microcrystals, and finite size ef- 
fects, so an extrapolation is used to establish a value for 
Bsat- Our ac susceptibility data shown in Fig. 5 suggest 
Bsat = 6.5 ± 0.2 T for T < 200 mK (when considering 
Bsat as the mid-point of the region where x{B) is most 
strongly changing), and this value agrees with the one 
extracted from the analysis of the Cm data at 9 T and 
14 T. 

The crossover region, ABgat) that spans from the re- 
gion from about 6 T to 7 T at low temperature, Fig. 5, 
merits further analysis. For example, when ksT « O.IJ, 
thermal smearing of the crossover is not expected. Al- 
ternatively, the powder character of the sample might in- 
duce a scatter of Bgat that can be evaluated as Bl^^—B^^^^ . 
Previous electron paramagnetic resonance studies pro- 
vided values for the anisotropic g-tensor of g|| = 2.25 and 
g± = 2.05.^^ Assuming an ideal square lattice and taking 
J/ks = 3 K, then 



-"sat 



4J 



5ii , Mb 



(10) 



yields | AB^at |< 0.8 T, which is less than but similar 
in magnitude to the observed width of the crossover. In 
addition, the finite size of the particles comprising the 
powder-like sample might cause broadening by limiting 
the in-plane correlation length ^ at the lowest tempera- 
tures. Perhaps more fundamentally, the spatial extent of 
^ near Bgat might be restricted by the underlying mag- 
netic frustration, which suppresses the antiferromagnetic 
correlations. 



C. Magnetic phase diagram of Cu(tn)Cl2 

Prior to building the magnetic phase diagram, the cor- 
relations between the shift of Tp and the position of high 
temperature specific heat maximum T^ax with respect 
to magnetic field are worth noting (Figs. 6 and 8). The 
monotonic increase of Tp with increasing magnetic field 
up to about 2 T corresponds to the relative field inde- 
pendence of Tmax in this field region. The relative field 



independence of Tp in the field from 2 to 4 T corresponds 
to the slight increase of Tmax, and finally, the rather rapid 
dccrciase of Tp above 4 T corresponds to a linear increase 

of Tj^ax- 

This relationship between the low temperature 
anomaly and the high temperature maximum suggests 
that the magnetic field mainly affects the short-range cor- 
relations. This interpretation is supported by the mono- 
tonic increase of Tp observed in low fields, and this be- 
havior is typical for the Berezinski-Kosterlitz-Thouless 
(BKT) transition theoretically predicted in the classical 
limit for low dimensional magnets in uniform magnetic 
fields.23'30 Quantum Monte Carlo studies of the 5=1/2 
isotropic square lattice in a magnetic field also revealed 
field-induced XY behavior, and a BKT transition at a 
finite temperature, Tbkt, was identified for sufficiently 
strong magnetic fields. Neither the presence of spatial 
anisotropy nor the introduction of the frustrated nnn 
coupling are expected to qualitatively change the physical 
picture derived for the isotropic square lattice. '^^ Conse- 
quently, the measured Tp vs. B dependence can be com- 
pared to the theoretical prediction of Tbkt vs. B calcu- 
lated for g = 2.12 and J/A'b = 3 K, and the comparison 
yields surprisingly good agreement between the theoreti- 
cal predictions and data (Fig. 8). It should be noted that 
the authors of Ref. 16 identified the BKT phase transi- 
tions at Tbkt, which lies below Tp. However, the deter- 
mination of Tbkt directly from the experiment is not as 
straightforward as the determination of Tp. From this 
point of view, one must be aware that the experimental 
determinations of the potential Tbkt are overestimated. 

Recent quantum Monte Carlo calculations of a i? vs. T 
magnetic phase diagram performed for the system of a 
tetragonal lattice with intralayer coupling J and inter- 
layer coupling J" revealed that, for sufficiently large spa- 
tial anisotropy J/ J" , the system preserves nonmonotonic 
B vs. T behavior typical for the BKT transition on the 
ideal square lattice. The enhancement of the critical 
temperatures and the saturation field with respect to the 
ideal 2D case is another effect of the intcrlaycr coupling. 
The predictions were compared with the experimentally 
established phase diagram of [Cu(HF2)(pyz)2]BF4, a rep- 
resentative of the quasi-2D Heisenberg square lattice with 
a finite temperature phase transition to 3D long-range 
order in zero magnetic field, and excellent agreement 
was found. -"^^ Unlike the phase diagram of Cu(tn)Cl2, 
the experimental data of [Cu(HF2)(pyz)2]BF4 arc sig- 
nificantly shifted above the theoretical prediction for the 
ideal square lattice model. The fact that our data lie 
significantly lower, i.e. below the theoretical predictions 
for the BKT transition (Fig. 8), supports the conjecture 
about the important presence of frustration. 

Considering all the ac susceptibility and specific heat 
features together, we can identify the major regions in 
the extended magnetic phase diagram for Cu(tn)Cl2 as 
a triangular magnet from the Neel phase (Fig. 9). At 
temperatures above 2 K, the system behaves as a para- 
magnet in all magnetic fields. The paramagnetic region 
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FIG. 9: (Color online) Extended magnetic phase diagram of 
Cu(tn)Cl2. The open squares represent the Tmax vs. B depen- 
dence, while the pluses (+) correspond to the line constructed 
from the condition that StS{T, B) - 5tS(T, 0) = 0, see text. 
The theoretically predicted BKT transition and the exper- 
imental Tp vs. B dependence are denoted by full triangles 
and full circles, respectively. Open circles correspond to the 
positions of the mid-points in the shoulders of ac susceptibil- 
ity data. The star represents the Bsat as extracted from an 
analysis of the specific heat data measured in 9 T and 14 T. 
The dotted and solid lines are guides for eyes. The regions 
of short-range order (SRO) and free or bound vortex (V) and 
antivortex (AV) pairs are identified. Detailed descriptions of 
each B — T region are given in the text. 



is separated by the line determined by the Tmax vs. B 
dependence. At temperatures below the line, 2D anti- 
ferromagnetic correlations develop, forming free vortices 
(V) and antivortices (AV) in the xy plane. The vor- 
tices are stabilized by the magnetic field, whose direc- 
tion defines the z-axis in the spin space. The forma- 
tion of bound V-AV pairs begins at temperatures below 
a line derived from the requirement that the difference of 
the entropy derivatives [StS{T, B) - StS(T, 0)] is zero.^^ 
This quantity equals the difference of the specific heat 
[C(r, B) - C{T, 0)] divided by temperature T. The in- 
tensity of the pairing process culminates at the temper- 
atures defined by the position of the low temperature 
specific heat anomaly. As was shown in the theoretical 
studies of the BKT transition, ^^'^^ the temperature of the 
BKT transition itself is about 20 — 30% lower than the 
position of the specific heat maximum. Consequently for 
Cu(tn)Cl2, the potential BKT transition can be expected 
at temperatures lower than the critical line constructed 
from the Tp vs. B data. 

Above the critical magnetic field region determined 
from the position of the high field shoulder in the isother- 
mal ac susceptibility data, out-of-plane ferromagnetic 
correlations induced by the magnetic field are stabi- 
lized, while in-plane spin correlations show paramagnetic 
behavior. '^^ In this region, a spin gap appears in the ex- 
citation spectrum, and this gap was detected by the ex- 



ponential character of the low temperature specific heat 
measured in 9 T. As shown in the phase diagram, the in- 
duced ferromagnetic correlations survive up to the tem- 
peratures on the critical line T^ax vs. B. Coincidence 
with the theoretical predictions for the square lattice, '^^ 
formation of the out-of-plane ferromagnetic correlations 
already begins at fields above 4 T as indicated by the 
linear dependence of Tmax vs. B above 4 T. 

The onset of quasi-long-range order (QLRO), induced 
by the field due to the stabilization of bound V-AV pairs 
expected for the unfrustrated square lattice at sufficiently 
low temperatures, should be reflected by the smaller en- 
tropy increase in finite magnetic field than in zero field. 
In other words, the presence of the V-AV pairs should 
lead to negative values of [C(T, B) — C{T, 0)] at low tem- 
peratures. However, as can be seen from the insets of 
Figs. 3 and 8, in Cu(tn)Cl2, the difference remains pos- 
itive down to the lowest temperatures. The increasing 
specific heat with magnetic field can be ascribed to the 
combined effect of the interlayer coupling and the exis- 
tence of soft spin modes, likely associated with the frus- 
trated nnn coupling, which can prevent the system from 
achieving full QLRO. 



V. SUMMARY 

The response of Cu(tn)Cl2 to an externally applied 
magnetic field has been investigated by specific heat 
and ac susceptibility studies at low temperatures, T > 
40 mK, and in magnetic fields up to 14 T. Specific fea- 
tures induced by the magnetic field in the behavior of 
both quantities allowed a magnetic phase diagram, whose 
regions are rooted in the presence of a high degree of 
short-range order, to be constructed. The quadratic tem- 
perature dependence of the specific heat observed below 
0.4 K and in magnetic fields up to, at least, 2.5 T arises 
from the predominantly 2D character of magnetic cor- 
relations and identifies a gapless excitation spectrum as 
predicted for the isotropic square lattice. Furthermore, in 
finite magnetic fields up to 6 T, a field induced anomaly 
appears near 0.8 K, and this feature is associated with 
the existence of a Berezinski-Kosterlitz-Thouless (BKT) 
transition theoretically predicted for the isotropic square 
lattice. The dimensionality of the observed magnetic 
field induced transition is unresolved and will be the fo- 
cus of future neutron and/or muon scattering studies. 
The absence of a phase transition in zero magnetic field 
is interpreted as a consequence of a combined effect of 
the weak interlayer coupling and the frustration within 
the magnetic layers. The saturation magnetic field value, 
as extracted from the specific heat and ac susceptibility 
data, deviates significantly from the one predicted for an 
isotropic square lattice, and this result is conjectured to 
be a consequence of the frustrating magnetic interactions, 
in the be layer, that have have been confirmed recently 
by quantum mechanical calculations.^'^ 

Finally, it is noteworthy that a dependence of the 
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specific heat has been reported for several frustrated 2D 
systems in zero magnetic field and was ascribed to various 
scenarios. While only a weak field dependence of the low 
temperature specific heat has been observed in Kagome^* 
and triangular^^ compounds, the linear increase of the b 
coefficient, as observed in Cu(tn)Cl2, points to the fact 
that even an infinitesimal field is able to introduce the 
changes. This sensitivity also supports the applicability 
of the BKT model to Cu(tn)Cl2; however, only to some 
extent. In the spin wave region, the magnetic field should 
decrease the specific heat and, correspondingly, the b co- 
efficient, as expected for the square lattice. Since the 
opposite tendency was observed, our experimental data 
suggest the presence of soft spin modes that are possibly 
connected with frustration and interlayer coupling. A 
spin wave analysis of the triangular magnet from the Neel 
phase in a magnetic field would be desirable to elucidate 
these observed differences. Recently, the BKT descrip- 
tion has been successfully used^^ to explain the 2D spin 
freezing transition observed in NiGa2S4, a model system 
for the S = I isotropic triangular Heisenberg lattice.'^'' 

In conclusion, the analysis of experimental data sug- 
gests that the 5=1/2 spatially anisotropic triangu- 
lar magnet from the coUinear Neel phase undergoes 



a Berezinski-Kosterlitz-Thouless transition induced by 
an applied magnetic field. Theoretical studies of the 
coUinear Neel phase in a magnetic field are necessary to 
specify the proper range of J, J' parameters where the 
BKT transition can occur. Microscopic magnetic studies 
of Cu(tn)Cl2, possibly employing neutron and/or muon 
scattering techniques, are needed to clarify the nature of 
the ground state in zero and nonzero magnetic fields. 
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